Unconventional resistivity at the border of metallic antiferromagnetism in N1S2 
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We report low-temperature and high-pressure measurements of the electrical resistivity p(T) of 
the antiferromagnetic compound MS2 in its high-pressure metallic state. The form of p(T) suggests 
that metallic antiferromagnetism in N1S2 is quenched at a critical pressure p c — 76 =b 5 kbar. Near 
p c the temperature variation of p(T) is similar to that observed in NiS2-ccSex near the critical com- 
position x — 1 where the Neel temperature vanishes at ambient pressure. In both cases p(T) varies 
approximately as T 1-5 over a wide range below 100 K. However, on closer analysis the resistivity 
exponent in N1S2 exhibits an undulating variation with temperature not seen in NiSSe (x = 1). This 
difference in behaviour may be due to the effects of spin-fluctuation scattering of charge carriers 
on cold and hot spots of the Fermi surface in the presence of quenched disorder, which is higher in 
NiSSe than in stoichiometric N1S2. 
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I. INTRODUCTION 

The electronic properties of metals on the border of 
magnetic phase transitions at low temperatures are of- 
ten found to exhibit temperature dependences that differ 
from the predictions of Fermi liquid theory. Early at- 
tempts to explore such non-Fermi liquid behaviour have 
been based on mean-field treatments of the effects of en- 
hanced magnetic fluctuations, as in the self-consistent 
renormalization (SCR) model. 1,2,3,4 

In a recent work the prediction of this model for the 
temperature dependence of the electrical resistivity p(T) 
was tested in a simple cubic d-metal, M3AI, at high pres- 
sures near to the critical pressure where ferromagnetism 
is suppressed^ The T 5 / 3 temperature dependence of the 
resistivity seen in M3AI and other related systems, where 
the magnetic correlation wavevector k(T) is small, ap- 
pears to be largely consistent with the SCR modelA In 
the idealized limit n — > at T — > 0, the SCR model 
predicts that in 3D the quasiparticle scattering rate 
varies linearly with the quasiparticle excitation energy, 
rather than quadratically as in the standard Fermi liq- 
uid picture. This form of is similar to that of the 
marginal Fermi liquid m.ode\^^ which is normally as- 
sociated with a linear temperature dependence of the re- 
sistivity. However, at the border of ferromagnetism the 
relevant fluctuations responsible for quasiparticle scat- 
tering are of long wavelength and thus are ineffective in 
reducing the current. This leads to a transport relax- 



ation rate r u 
but as T 5 / 3 



that differs from r q ^ 



and varies not as T, 



In practice, this form of r^ 1 or of the resistivity p(T) 
can be restricted to a relatively narrow range of temper- 
atures and pressures where k(T) is small compared with 
the characteristic wavevector of thermally excited mag- 
netic fluctuations. The SCR model is not restricted to 
this limit and can in principle include the effects of k(T) 
and even determine k(T) in a self-consistent fashion. We 
note that the SCR model assumes implicitly that an effec- 
tive underlying mechanism exists to remove momentum 
from the electron system. The validity of this assumption 
has not been clearly confirmed theoretically, but seems 
to be consistent with experiment in the cases mentioned 
above. In its conventional form the model does not in- 
clude the possible effects of inhomogeneities or texture 
that may arise on the border of first order ferromagnetic 
transitions as in, e.g., MnSi^ii 4 ^ The SCR model is 
also expected to fail on the border of electron localiza- 
tion, as near a Mott transition or close to local quantum 
critical points in heavy electron compounds™^ 

The applicability of the SCR model has also been ques- 
tioned for the case of itinerant-electron antiferromag- 
netism in general, even well away from the border of a 
Mott transition^ii2i2S In this paper we present an at- 
tempt to test the prediction of the SCR model in a metal 
on the border of metallic antiferromagnetism in 3D for 
which p is predicted to vary as T 3 / 2 in the idealized limit 
n — > 0, T — > 0, where n now stands for the correlation 
wavevector for the staggered magnetization. The expo- 
nent of 3/2 is the ratio of the spatial dimension d = 3 
and the dynamical exponent z = 2. This may be con- 
trasted with the corresponding exponent of 5/3, which 
is the ratio of (d + 2) = 5 and z = 3, for the border 
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FIG. 1: Temperature- composition phase diagram of 
NiS 2 -x Se^."*" Lower inset shows the temperature depen- 
dence of the resistivity of NiSSe (x — 1), which is just 
on the border of a metallic antiferromagnetic state at low 
temperatures J- The upper inset is the pyrite crystal struc- 
ture of NiS2 and NiSe2. PI and PM stand for paramagnetic 
insulator and paramagnetic metal, respectively. AFI and 
AFM stand for antiferromagnetic insulator and antiferromag- 
netic metal, respectively. SR-AFI stands for short range an- 
tiferromagnetic insulator The temperature-pressure phase 
diagram of MS2 is expected to be similar in form to the 
temperature-composition phase diagram of NiS2-xSe cc . In 
this study, the AFM boundary of MS2 is found to be at 
p c = 76 =b 5 kbar. Thus, NiS 2 at p c and NiSSe at ambient 
pressure may be expected to be similar except for the level of 
quenched disorder, which is normally higher in NiSSe than in 
NiS2- The residual resistivity of NiSSe is typically an order 
of magnitude higher than in MS2. 



of metallic ferromagnetism. (The 2 in d + 2 arises from 
the effect of small-angle scattering that is absent in the 
case of the staggered magnetization.) This simple model 
for the scattering from antiferromagnetic fluctuations as- 
sumes that the scattering rate can be averaged over the 
Fermi surface. Within the SCR model this procedure 
would seem to require the presence of a sufficient level 
of quenched disorder assumed to have only the simple 
consequence of inhibiting the short circuiting caused by 
the carrier on the cold spots of the Fermi surface, i.e., 
regions far from the hot spots connected by the antiferro- 
magnetic ordering wavevector and thus strongly affected 
by spin-fluctuation scattering 

The effect of quenched disorder on the temperature de- 
pendence of p in the SCR model shows up particularly 
clearly in the temperature-dependent resistivity expo- 
nent denned as n = d In Ap/d In T, where Ap = p(T) — po 
and po is the residual resistivity, i.e., the resistivity ex- 
trapolated to T = K. The resistivity exponent n may 
be described in terms of the Fermi wavevector hp, the 
elastic mean free path / of charge carriers and the re- 
duced temperature t = T/T s f, where T s f is a character- 
istic spin fluctuation temperature^ With increasing t, n 
drops from 3/2 towards unity around t w l/(fcpZ), back 



towards a value of the order of 3/2 around t ~ 1/y/kpl 
and then towards zero for t » l/(fe F /) j 18 i 19 i 2Q This 
undulating behaviour of n is a dramatic prediction of 
the model for relatively pure samples that could in prin- 
ciple be tested by studying a series of samples of the 
same material with different levels of quenched disorder. 
RoschiSiSl introduced this model in an effort to under- 
stand the temperature dependence of the resistivity ex- 
ponent measured for the f-electron metal CePd2Si2^i*2^ 
The results were, in this case, not entirely conclusive be- 
cause of the anisotropic and complex nature of the spin 
fluctuation spectrum in this material. 

Here we present a test of Rosch's model in a d-electron 
system with a cubic structure and T s f one to two orders 
of magnitude greater than in typical heavy f-electron sys- 
tems. We compare the temperature dependence of p of 
two closely related systems on the border of antiferro- 
magnetism. The two materials are NiS 2 near 76 kbar and 
NiS2-xSe^ for x = 1 where the Neel temperature T/v van- 
ishes at ambient pressure. Due to random variations in 
the locations of the S and Se atoms in the lattice, values 
of po found in NiSSe are typically one order of magnitude 
higher than in stoichiometric MS2 compounds. 

MS2 crystallizes in the cubic pyrite structure and 
is an antiferromagnetic insulator at ambient pres- 
sure at low temperatures ^iS4 It can be metallized 
via the application of pressur o 25 i 26 i 27 i 28 or by Se 
substitutioni 10 i 12 i 25 i 27 i 28 i 29 i 30 The temperature-pressure 
phase diagram of MS2 and the temperature-composition 
phase diagram of NiS2-^Se ;r (Fig. [IJ are expected to 
be similar. Of particular interest for this paper is the 
boundary of metallic antiferromagnetism that appears 
(i) at x = 1, i.e., for NiSSe, in temperature-composition 
phase diagram (Fig. and (ii) at p c in the temperature- 
pressure phase diagram for stoichiometric MS2. We note 
that in both cases the quantum critical point for metallic 
antiferromagnetism arises well beyond the metal insula- 
tor transition (see Fig. for the case where composition 
is the control parameter). 

We present a high-pressure study of the temperature 
dependence of the resistivity of MS2 in the metallic state, 
which reveals a critical pressure p c = 76 kbar. The results 
at p c are compared with that of NiSSe at ambient pres- 
sure reported in Ref. 0. In both materials one observes 
a non- Fermi liquid form of the resistivity that appears in 
first approximation to be consistent with the predictions 
of the SRC model (inset of Fig. Q for NiSSe). However, 
the resistivity of MS2 near p c has an undulating compo- 
nent in its variation with temperature that is not seen in 
NiSSe. This difference in behaviour might arise from the 
effects of cold and hot spots on the Fermi surface as sug- 
gested by the model due to Rosch discussed above. 19,20 



II. EXPERIMENTAL 

Single crystals of MS2 have been grown by chem- 
ical vapour transport using iodine as the transport 
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agent. 25 The residual resistivity ratio (p(273 K)/po) 
above 50 kbar is about 30 for our samples, compared 
with about 3 for NiSSe, 12 The carrier mean free path of 
MS2 is thus expected to be about an order of magnitude 
higher than in NiSSe. 

Pressure was applied by means of a non-magnetic 
Bridgman cell using tungsten carbide anvils. The cell 
used is a scaled down version of that designed by 
Wittig^i The cullet diameter of the anvils was 3.5 mm. 
The gasket was made of pyrophyllite with a central hole 
to accommodate the NiS 2 sample as well as a Pb sam- 
ple and the steatite pressure-transmission medium. Af- 
ter compression, the sample space was reduced to about 
1.5 mm in diameter and 0.1 mm in thickness. The pres- 
sure was determined from the superconducting transition 
temperature of the Pb sampled The width of the super- 
conducting transition suggested that the pressure varia- 
tion over the Pb sample was about 10% of the average 
pressure. 

The resistivity was measured via the 4-terminal ac 
technique. Four 50 pm Pt leads were passed into the 
high-pressure region via grooves in the insulating pyro- 
phyllite gasket. The bare wires were rested on top of 
the samples and pressed onto the sample surface during 
pressurization to achieve adequate electrical contacts. 

The resistivity measurements were carried out with 
two different and independent systems, a helium circula- 
tion cryostat (ILL Orange cryostat; 1.5 K to 300 K) in 
Grenoble and an adiabatic demagnetization refrigerator 
(0.04 K to 100 K) in Cambridge. The latter system had 
two voltage channels, one with a low temperature trans- 
former and the other with a room-temperature trans- 
former. The excitation currents were 1 mA and 160 pA 
in the Orange cryostat and in the adiabatic demagnetiza- 
tion fridge, respectively. The results obtained with these 
two experimental systems were consistent with each other 
where comparisons could be made. 



III. RESULTS 

The temperature variations of the resistivity of MS2 
in the high pressure metallic state above 40 kbar are pre- 
sented in Figures El and 03 (The insulator-to- metal tran- 
sition (not shown) was observed at around 30 kbar, in 
agreement with the literature^ 5 -) Figure shows Ap vs 
T up to 80 K and Figure 03 is an expanded view of Ap vs 
T in the range 0.05 K to 2 K. Plots of Ap vs T 2 and Ap 
vs T 3 / 2 over different ranges in temperatures are shown, 
respectively, in Figures [5] and [SI 

Below 1 K the resistivity can be described by an equa- 
tion of the form p = p + AT 2 over the entire pressure 
range explored, 43 kbar to 96 kbar. The pressure varia- 
tions of the fitted values of A and po are given in Figure0] 
The T 2 coefficient A exhibits a peak at p c = 76 ± 5 kbar. 

Figures and El compare Ap vs T 2 and Ap vs T 3 / 2 , re- 
spectively, in three panels each covering different ranges 
in temperature. Figure [St, in particular, highlights the 




FIG. 2: Temperature dependence of the resistivity in the 
metallic state of MS2 at high pressures. Ap is p(T) — po, 
where po is the residual resistivity extrapolated to T — K. 
The curves are shifted vertically for clarity. 



0.1 



0.075 



£o.05 

Q_ 

< 

0.025 



(a) 

. p = 76 kbar 
72 kbar 
66 kbar 
63 kbar 
57 kbar 
43 kbar 

A P = P"Pn . 



NiS 



76 kbar #/' 



43 kbar 




0.1 



0.075 



0.05 



Q_ 

< 



0.025 



(b) 

p = 76 kbar 
79 kbar 
83 kbar 
88 kbar 
91 kbar 
96 kbar 



NiS 



76 kbar 



s 




96 kbar 



1.5 



FIG. 3: Low temperature variation of the resistivity in the 
metallic state of NiS 2 at high pressures. Ap grows in strength 
up to 76 kbar (a) and weakens gradually above 76 kbar (b). 



0.03 



| 0.02 

o 

a 

zL 



0.01 



NiS 




a 
it 



P = P + AT 



30 p(kbar) 110 
fit from 50 mK to 1 K 



0L 
30 



50 



70 
p (kbar) 



90 



110 



FIG. 4: The T 2 coefficient of the resistivity, A, and residual 
resistivity, po (inset), of MS2 in the metallic high-pressure 
state. The parameters A and po are obtained by a fit of 
the resistivity as indicated in the figure. A is peaked at 
p c = 76 db 5 kbar. This marks the boundary of the metal- 
lic antiferromagnetic state of NiS2- 



T 2 variation observed at all pressures with a peak of A at 
around p c as discussed above. Figures Et and0 on 
the other hand, suggest that near this pressure Ap varies 
roughly as T 3 / 2 over a decade in temperature above a few 
K. This is the behaviour predicted by the SCR model for 
a system on the border of metallic antiferromagnetism 
as discussed in the introduction. The existence of a T 2 
regime even at p c , however, suggests that the transition 
into the antiferromagnetic state below p c may not be con- 
tinuous, i.e., that the antiferromagnetic quantum critical 
point is not quite reached due to the onset of a first order 
transition. We note that the weak pressure variation of 
the Fermi liquid crossover temperature Tfl near p c (see 
Fig. and the caption) would seem to rule out an expla- 
nation of the T 2 resistivity in terms of an inhomogeneity 
in pressure. 

The identification of p c with the boundary of metal- 
lic antiferromagnetism is also suggested by the corre- 
spondence of the temperature-pressure and temperature- 
composition phase diagram of MS2 found in earlier 
wor k Jfiil2i22i2i This predicted that the critical pressure 
for the border of metallic antiferromagnetism should be 
around 60 kbar, which is of the same order of magni- 
tude as p c denned above. The identification of p c with 
the antiferromagnetic boundary in the metallic state is 
tentative and needs confirmation by other measurements 
and in particular the detection of a signature of T/v in the 
resistivity for pressures below p c . Here we focus atten- 
tion mainly on the behaviour of p near p c and contrast it 
with that of NiSSe at ambient pressure (inset of Fig.^J. 

The non- Fermi liquid behaviour of Ap over a wide 
temperature range near p c (low-temperature data at 
76 kbar with added 77 kbar data up to 300 K; the high- 
temperature data is scaled to match the low-temperature 
data) is highlighted in Figure [SI Panel (a) shows the 
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FIG. 5: Resistivity vs T 2 in the high pressure metallic regime 
of NiS 2- The three panels are for different temperature in- 
tervals. A quadratic temperature dependence of Ap is seen 
only at the lowest temperatures (below a few K) and the T 2 
coefficient of Ap is peaked at p c ~ 76 kbar as shown in Fig.0] 



dramatic upturn of Ap/T 2 , which in the Fermi liquid 
regime is expected to saturate to a constant value. Panel 
(b) is a plot of In Ap vs InT which shows that the aver- 
age slope corresponds to a resistivity exponent close to 
3/2, as discussed above and as seen in NiSSe at ambient 
pressure. However, in contrast to NiSSe, the tempera- 
ture dependence of the resistivity in MS2 at p c exhibits 
an undulating structure which is evident in Figure |S| and 
highlighted in the temperature dependence of the resis- 
tivity exponent n = SlnAp/<91nT shown in the inset. 
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FIG. 6: Resistivity vs T 3 ^ 2 in the high-pressure metallic 
regime of MS2. The three panels are for different temper- 
ature intervals. An approximately T 3 ^ 2 variation of Ap is 
seen around p c ~ 76 kbar only over a decade in temperature 
above a few K. Ap is quadratic in temperature below a few K 
at all pressures in the metallic regime studied (see Fig- Et) • 

We note that this undulating structure has been observed 
in several samples of N1S2 and in two independent mea- 
surement systems. 

IV. DISCUSSION 

The resistivity measurements suggest that antiferro- 
magnetism in the high-pressure metallic state of N1S2 is 
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FIG. 7: Contours of the resistivity exponent n = 
d\n Ap/d\n.T in the temperature-pressure plane for N1S2 in 
the metallic regime. A Fermi liquid (FL) T 2 temperature de- 
pendence of the resistivity is seen at the lowest temperatures 
up to Tfl defined, somewhat arbitrarily, by the condition 
ti(Tfl) = 1.8. Tfl reaches a minimum at p c ~ 76 kbar. A 
non- Fermi liquid T n form of the resistivity with n ~ 3/2 is 
seen over approximately a decade in temperature above a few 
K near p c . 



suppressed at a critical prssure of p c = 76 ± 5 kbar. At 
p Cl we find that the T 2 coefficient of the resistivity A has 
its maximum and the Fermi liquid crossover temperature 
Tfl defined in the caption of Figure Q has its minimum. 
Above Tfl near p c the resistivity has an approximately 
T 3 / 2 temperature dependence over a decade in temper- 
ature. As already stated, this is the behaviour expected 
in the SCR model for a metal in 3D on the border of 
antiferromagnetism at low temperatures. Also, we note 
that p c is of the same order of magnitude as that in- 
ferred for the correspondence between the temperature- 
composition and temperature-pressure phase diagram as 
discussed in the previous section. 

The weakness or absence of a signature of Tat in the 
resistivity below p c is not necessarily surprising since the 
same is the case for NiS2- ;E Se :E in the range 0.4 < x < 1 
where antiferromagnetic order has been con- 
firmed by neutron scattering iS^^ The signature of T/y 
in p can be washed out by inhomogeneities in composi- 
tion or pressure and may only show up in high-precision 
measurements of dp/dT in the pressure range where Tn 
is not too strongly varying with pressure and where pres- 
sure is hydrostatic. 

The fact that Tfl and correspondingly the T 2 coeffi- 
cient of the resistivity A remain finite at p c is not nec- 
essarily inconsistent with our assumption that p c marks 
the border of antiferromagnetism. The antiferromagnetic 
transition vs pressure may be first order as is found in a 
number of antiferromagnetic metals such as YMn2 and 
GdMn2™ First order quantum phase transitions also 
seem to be common in ferromagnetic metals such as MnSi 
and Ni 3 Ah 5 i 37 i 38 T F l and A remain finite at the anti- 
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FIG. 8: The non- Fermi liquid form of the resistivity of NiS2 
near p c « 76 kbar. (a) Ap/T 2 vs T does not saturate and 
thus has a non-Fermi liquid form over a wide temperature 
range, except at very low temperatures (Fig.0. (b) A plot of 
In Ap vs In T and the resistivity exponent n = d In Ap/d In T) 
vs InT exhibits an undulating form in contrast to the sim- 
ple T 3/2 temperature dependence seen in NiSSe (Ref.IH and 
Fig-E- This undulating form of the resistivity has been seen 
in different N1S2 specimens and with two independent mea- 
surement systems. 



ferromagnetic boundary of YMn2 , a system that shares 
some features in common with MS2. (We note, how- 
ever, that the absence of Fermi liquid behaviour does not 
necessarily mean that the quantum phase transition is 
continuous 

The Fermi liquid crossover temperature, TpLi m NIS2 
is very much lower than the characteristic spin fluctu- 
ation temperature, T s /, which may be expected to be 
of the order of 10 3 K in a typical d-metal. We con- 
sider to what extent the behaviour above Tfl, but be- 
low T s j, may be understood in terms of an itinerant- 
electron model for a continuous antiferromagnetic quan- 
tum critical point. The predictions of the SCR model 
as analyzed by Rosch has already been outlined in the 
introduction for an antiferromagnetic quantum critical 
point. With increasing reduced temperature t = T/T s f, 



the SCR model predicts that the resistivity exponent 
n = <91nAp/<91nT drops from an initial value of 3/2 
towards unity around t ~ l/(fcpZ), back to a value of 
order 3/2 around t « and then to zero for 

t » l/(fc j pZ) JSil2i2Si This behaviour is qualitatively sim- 
ilar to that seen in MS2 at p c and above Tfl (inset 
Fig. EJ). For reasonable choices of parameters for MS2, 
T 3f « 10 3 K, k F « 0.5 A" 1 and I « 140 A" 1 , 39 we expect 
the minimum of n to occur near 15 K and the maximum 
at around 120 K. These crossover temperatures are in 
rough agreement with our observations. 

In NiSSe, I is an order of magnitude smaller than in 
MS2 and thus the minimum and maximum of n are ex- 
pected to arise at around 100 K and 320 K, respectively. 
Over the temperature range of the experiments shown in 
the inset of Figure ^ the model predicts a simple T 3 / 2 
temperature dependence without modulation, as is seen. 
The difference in the behaviour of N1S2 and NiSSe at 
their respective critical conditions is thus not surprising. 
The combined effects of scattering from spin and lattice 
fluctuations might also lead to an undulating behaviour 
of n, but this would naively be expected to arise in both 
MS2 and NiSSe, in contradiction with observation. 



V. CONCLUSIONS 

The temperature dependence of the resistivity of MS2 
near the critical pressure, which has been found to be 
p c = 76 ± 5 kbar, is consistent with that expected 
for a metal on the border of itinerant-electron antifer- 
romagnetism. The Fermi liquid crossover temperature, 
Tfl, which defines the range over which the resistivity is 
roughly quadratic in temperature, does not vanish at p c , 
but is three orders of magnitude smaller than the charac- 
teristic spin fluctuation temperature T s f . The finite value 
of Tp l may indicate that the antiferromagnetic quantum 
critical point is first order, as is the case in related ma- 
terials such as YMn2. 

Over a wide temperature range above Tfl the tem- 
perature variation of the resistivity exhibits a non-Fermi 
liquid form that can be understood in terms of the effects 
of spin-fluctuation scattering at cold and hot spots of the 
Fermi surface as anticipated by Rosch in his refined treat- 
ment of the self-consistent-renormalization (SCR) model. 
In particular, the resistivity exponent n = din Ap/dlnT 
exhibits an undulating structure which is consistent not 
only qualitatively, but also approximately quantitatively 
with the predictions of this model. The model also ac- 
counts for the absence of this undulating structure in the 
related, but less pure material, NiSSe, which is known to 
be at the border of antiferromagnetism at ambient pres- 
sure. 
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